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Observational evidence for the link between the variable 
optical continuum and the subparsec-scale jet of the radio 
galaxy 3C 390.3 
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ABSTRACT 

The mechanism and the region of generation of variable continuum emission are poorly 
understood for radio-loud AGN because of a complexity of the nuclear region. High- 
resolution radio VLBI (very long baseline interferometry) observations allow zooming 
into a subparsec-scale region of the jet in the radio-loud galaxy 3C 390.3. We com- 
bined the radio VLBI and the optical data covering the time period of 14 years to 
look for a link between optical flares and parsec-scale jet. We identify two stationary 
and nine moving radio features in the innermost subparsec-scale region of the jet. All 
nine ejections are associated with optical flares. We found a significant correlation 
(at a confidence level of > 99.99%) between the ejected jet components and optical 
continuum flares. The epochs at which the moving knots pass through the location 
of a stationary radio feature and the optical light curve reaches the maximum are 
correlated. The radio events follow the maxima of optical flares with the mean time 
delay of 0.10 ± 0.04 years. This correlation can be understood if the variable optical 
continuum emission is generated in the innermost part of the jet. A possible mech- 
anism of the energy release is the ejection of knots of high-energy electrons that are 
accelerated in the jet flow and generate flares of synchrotron continuum emission in 
the wide range of frequencies from radio to X-ray bands. In this scenario, the beamed 
optical continuum emission from the jet and counterjet ionizes a gas in a subrelativis- 
tic outflow surrounding the jet, which results in a formation of two outflowing conical 
regions with broad emission lines (in addition to the conventional broad line region 
around the central nucleus). 

Key words: galaxies: jets - galaxies: nuclei - galaxies: individual: 3C 390.3 - radiation 
mechanisms: non-thermal. 



1 INTRODUCTION 

The variable continuum flux in AGN, signaling the activ- 
ity of the central engine, is detected throughout the entire 
electromagnetic spect r um, on time-scales from days to years 
dPeterson etHdl |2002|; IZhend Il996l ; IWamsteker et all Il997l : 
IShapovalova et ail l200ll ). The bulk of the continuum flux 
is believed to be generated in the accretion disk and is re- 
sponsible for ionizing the cloud material in the broad-line 
region (BLR). Localization of the source of the variable con- 
tinuum emission in AGN is instrumental for understanding 
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the mechanism for release and transport of energy in active 
galaxies. In radio-quiet AGN, representing about 90 % of 
the AGN population, the presence of rapid X-ray flux vari- 
ations and iron emission line (Fe Ka) indicates that most 
of the soft X-ray emissio n originates from the accretion disk 
l|Mushotzkv et alj Il993h . In radio-loud AGN, the activity 
of the central engine is accompanied by highly-relativistic 
collimated outflows (jets) of plasma materi al formed an d 
accelerated in the vicinity of the black hole (|Ferraril fl998'). 
Inhomogeneities in the jet plasma appear as a series of com- 
pact radio knots (jet components) observed on s cales rang- 
ing from several light weeks to about a kiloparsec (|Alef et al.l 
Il99rj ; iKellermann et al.ll2004f ). 
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Figure 1. Radio structure of 3C 390.3 observed in 1998.21 with 
very long baseline interferometry at 15 GHz (2 cm). Innermost 
fraction of the jet is shown in the inset. The resolving point-spread 
function (beam) plotted in the upper left corner is 0.87 mas X 0.55 
mas oriented at an angle of 8° (clockwise rotation). The peak 
flux density in the image is 190mJy/beam (3.1 X 10 9 K) and the 
rms noise is 0.2 mjy/beam. The contours are drawn at 1, V2, 2 ... 
of the lowest contour shown at 0.6 mjy/beam. The labels mark 
three stationary features (D, SI and S2) and a subset of moving 
components (C2— C6) identified in the jet. Note that six more 
components, C7-C12, have been first identified in the jet in the 
VLBA images at later epochs (see Fig. [2}. 



The unification scheme (|Urrv fe Padovam] Il995t ) of 
radio-loud AGN suggests that the optical continuum as- 
sociated with the accretion flows and broad emission lines 
are viewed directly in radio-loud quasars and BL Lacertae 
objects (BL Lacs). The Doppler boosted continuum emis- 
sion from the re lativistic jet may dominate at a ll energies 
in BL Lacs (see lUlrich et al.lll997l ; IWorralll l2005t ). suggest- 
ing that the continuum variability in radio-loud AGN is re- 
lated to both the jet and the instabilities of acc retion-disk 
flows (|Mushotzkv et al.lll993l ; lUlrich et ai1ll997f ). In radio- 
loud quasars, the beaming of continuum emission is less 
but opti cal synchrotron r adiation from the core of the jet 
(3C 273; ISoldi et al.ll200sl ) can be still significant. In radio- 
loud galaxies the emission from accretion disk or hot corona 
can be hidden by an obscuring torus, and hence, the bulk 
of the optical continuum and broad-line emission in radio 
galaxies can be attributed to the relativistic jet rather than 
the central engine if the boosted jet-emission is still strong, 
i.e. the jet is intrinsically luminous and highly relativistic. 

The presence of a positive correlation between beamed 
radio luminosity of the jet and optical nuclear luminosity in 
the sample of radio galaxies suggests that the optical emis- 
sion is non-thermal and may o riginate from a relativistic je t 
l|Chiaberge et ah 1 1 19991 . 2002; H ardcastle fc Worralll [2000 ). 
In radio-loud quasars, there is evidence for a variable syn- 
chrotron optic al flare generated around the radio core of 
PKS 0420-014 (|D'Arcangelo et al.l I2007T ) an d in the inner- 
most ~ 0.2 pc region of the jet in BL Lac (Marschcr et al. 
2008). A correlation between ejection of superluminal knots 
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Figure 2. Separation of the jet components relative to the sta- 
tionary feature D (blue circles) for 21 epochs of VLBI observa- 
tions. The stationary components SI and S2 are marked by open 
red circles. Designations of moving components are shown in the 
upper left corner. The lines represent the best linear least-squares 
fits to the component separations. 



and dips in the X-ray emission was re ported for the ra- 
dio galaxy 3C 120 (jMarscher et al.ll2002h . Evidence for cor- 
related radio and optical variations has been reported for 
several AGN with time lags up to few hundreds days 
ijHanski et al.l l2002). The detection of a correlation between 
variability of the continuum flux and changes in the ra- 
dio structure in a radio galaxy on subparsec-scales would 
be the most direct evidence of optical continuum emission 
coming from the jet. We combine the results from moni- 
toring of the double-peaked broad line radio-loud galaxy 
3C 390.3 (redshift z = 0.0561) in the optical from 1992 
to 20 08 time period (|Shapovalova et al1l200ll ; ISergeev et alj 
2002, Shap ovalova et al, in preparation), UV |Zhendll996 
and X-ray l)Leighlv et al ] 1 19971 : iGliozzi et alJl2oW regimes 
with 21 VLBI observations of its radio emission made 
from 1994 to 2008 as p art of the 15 GHz VLBA0 survey 
ijKellermann et al.l 120041 ) and our VLBA observations dur- 
ing the time period from 2005 to 2008. 

In Sect. [2] we analyze the structure and kinematics of 
the pc-scale jet in 3C 390.3. Correlations between properties 
of the compact jet and nuclear optical emission on scales 
less than one parsec are analyzed in Sect. [3] Variable ra- 
dio emission from stationary components of the jet is dis- 
cussed in Sect. [4] and the identification of these stationary 
components is discussed in Sect. [S] In Sect. [S] we discuss 
possible radiation mechanism acting in the nuclear region of 
3C 390.3. 



2 STRUCTURE AND KINEMATICS OF THE 
PARSEC-SCALE JET 

To parameterize the structure of the radio emission , we ap- 
plied the technique of modelfitting (|PearsorJll997r i and fit 
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Table 1. Kinematics and ejection epochs of moving jet components in 3C 390.3. 



Comp. 


to 


tsi 




(yr) 


(yr) 


(1) 


(2) 


(3) 


C4 


1992.79±°,'g 


1993.70l°'3g 


C5 


1995.62+.°;°^ 


1996.21+.° : ° 4 


C6 


1995.66t°22 


1996.75+.°;°? 


C7 


1997.15±g:|f 


1998.43+.°;^ 


C8 


1999.29+°;^ 


2000.38+°;}° 


C9 


2001.651°;" 


2002. 71+. 


CIO 


2003.11±° : || 


2004.17t°'29 


Cll 


2004.89^°;^ 


2005.84l°j 2 


C12 


2006.631°;^ 


2007.73l°'g4 



Mr 

(mas yr -1 ) 
(4) 



/^app 

(5) 



0.29±0.03 
0.46±0.03 
0.24±0.02 
0.21±0.03 
0.22±0.04 
0.26±0.05 
0.25±0.04 
0.28±0.06 
0.28±0.08 



1.05±0.11 
1.63±0.11 
0.87±0.09 
0.74±0.09 
0.78±0.13 
0.90±0.16 
0.90±0.14 
1.01±0.21 
1.00±0.27 



Columns are as follows: (1) - Component identifier, (2) - Ejection time from the component D, (3) - 
Passing time through the stationary component SI, (4) - Angular radial speed, (5) - Radial speed in units 
of the speed of the light. 



interferometric visibility amplitudes and phases from each of 
the 21 VLBA datasets by a set of two-dimensional, circular 
Gaussian features (shaded circles in Fig. [I}. Similar fits have 
been obtained for 17 observations from t he 15 GHz VLBA 
survey database l|Kellermann et al"1l2004l ) and for four our 
observations with the VLBA at 15 GHz, for the purpose of 
cross-identifying and tracing different features in the jet (Ta- 
ble [Al] in the Appendix A). We use the positions and flux 
densities of these components for tracing the evolution of 
the jet emission on angular scales of ~ 3 milliarcseconds 
(mas). At the distance of 3C 390.3, lmas corresponds to a 
linear distance of 1.09 pc for the flat ACDM cosmology with 
the Hubble constant Ho = 70 km s _1 Mpc -1 and the matter 
density fi m = 0.3. 

Based on the separations, sizes and flux densities of 
the Gaussian components located within 3 mas of the fea- 
ture D at the narrow end of the jet (Fig. [T]), we identified 
nine moving components C4-C12 (in addition to the pre- 
viously known components C2 and C3; Alef et al. 1996) 
and two stationary features SI and S2 separated from D by 
(0.28±0.03) mas and (1.50±0.12) mas, respectively (Fig.0. 
Note that two innermost components D and SI are spatially 
resolved and their flux densities are uncorrelated (see Ap- 
pendix A). The more distant stationary feature S2 is prob- 
ably related to a small change in the direction of the flow 
that causes the re l ativis tic brightening of the radio emis- 
sion (|G6mez et al.lll997i ). The presence of SI may be re- 
lated to both geometrical and physical factors. Recent re- 
sults irorri_£k2s^_yLBImonitoring of a sample of parsec-scale 
jets (jjorstad et al.|[200lf ) imply that such stationary features 
are likely to be common in compact relativistic flows. 

Proper motions of (0.2 to 0.5) mas/yr (which corre- 
spond to apparent speeds of (0.7 to 1.6) c, where c is the 
speed of light) are estimated from linear fits to the observed 
separations of the components C4-C12 from the component 
D (Table [TJ. Resolution of the VLBA data at 15 GHz is 
not sufficient for detecting moving components at locations 
between D and SI, since D and SI are strong, unresolved 
features separated by less than one beamwidth, while the 
moving features are less prominent and more extended. High 
frequency VLBI observations are required to study these an- 



gular scales in 3C 390.3. Back-extrapolation of the fits for C5 
and C6 indicates that these two components may originate 
from a single event. They could result from a moving per- 
turbation i n the jet that crea tes a forward and a reverse 
shock pair ijGomez et al.lll997|) or trailing features behind 
a stro ng relativistic shock ( Gomez et al.|[200l] ; lAgudo et al.l 
l200ll ). We use the linear fits to estimate, for each moving 
component, the epoch, to, at which it was ejected from the 
component D and the epoch, tsi, when it passed through the 
location of the stationary feature SI (see Tabled] Figs.[2]and 
EjJ. For C12 component, these epochs are estimated using the 
mean slope of the best-fits of eight components (C4-C11). 



3 LINK BETWEEN SUBPC-SCALE JET AND 
VARIABLE OPTICAL CONTINUUM 
EMISSION 

Optical fluxes and its errors measured at 5100A (rest frame) 
from observations of 3C 390.3 between 1995 to 1 9 99 (F ig. |3] 
red squares) are pub lished in IShapovalova et all l|200ll l and 
ISergeev et ail <|2002l ). During the time period from 2000 to 
2007, spectra of 3C 390.3 (97 nights) were taken with the 
6 m and 1 m telescopes of the Special Astrophysical Obser- 
vatory of Russian Academy of Sciences (SAO RAS, Russia) 
and at the 2.1 m telescope of the "Guillermo Haro Observa- 
tory" at Cananea, Sonora (Mexico), respectively. Long-slit 
spectrographs, equipped with CCD detector arrays, were 
used to cover the typical wavelength interval from 4000 A 
to 7500 A with spectral resolution between 5 A and 15 A 
and the signal-to-noise ratio >50 in the continuum near Ha 
and H/3 emission lines. The uncertainty in our flux deter- 
minations for the continuum is <3 %. Detailed description 
of the method and measure ments of optical fluxes are given 
in lShapovalova et al.l l|200ll ). Optical continuum fluxes mea- 
sured at 5100 A between 2000 and 2007 are presented in the 
Table [2] (for a full table see the online Supporting Informa- 
tion). 

Variations of the optical continuum emission (Fig. [3] red 
squares) exhibit a superposition of weak flares (few weeks 
to several months) , prominent flares (from months to years) 
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Figure 3. The varia tions of the continuum fluxes of 3C 390.3 du ring the 1992 to 2008 ti me period. O ptical continuum light curve at 
5100A (red squares; Sh apovalova et alj feoOl: Scrgc evet al.l|20 02j). UV continuum fluxes |zhcng 1996) scaled by a factor of 50 (grey 
squares), and variations of radio flux density from D and SI components (filled cyan diamonds and circles) are presented. The times of 
ejection, trj, of radio knots from D and the times of their separation, tsii from the component SI are marked by open triangles and filled 
triangles respectively. 



Table 2. The first ten entities of optical continuum fluxes of 
3C 390.3 from the online electronic table. Columns: 1 - UT date; 
2 - Julian date; 3 — continuum flux at 5100 A (in units of 10~ 15 
erg s — 1 cm _2 A _1 ). 



UT-date JD Flux at 5100 A 

(2400000+) 



1 


2 








3 




2000Jull9 


51745 


363 


1 


557 


± 


0.047 


2000Jul22 


51748 


230 


1 


602 


± 


0.048 


2000Jul30 


51756 


293 


1 


543 


± 


0.046 


2000Oct05 


51823 


141 


1 


410 


± 


0.042 


2000Novl8 


51867 


137 


1 


558 


± 


0.047 


2000Novl8 


51867 


223 


1 


661 


± 


0.050 


2000Novl9 


51868 


129 


1 


660 


± 


0.050 


2001Janl6 


51925 


633 


1 


533 


± 


0.046 


2001Jan20 


51929 


621 


1 


465 


± 


0.044 


2001Marl3 


51981 


949 


1 


396 


± 


0.042 



and a long-term rising trend (few decades). The later is seen 
also in the histori cal light curve of the op tical continuum 
in 3C 390.3 (e.g. IShapovalova et all 12001). It is remark- 



able, that all nine ejections are associated with nine opti- 
cal flares happening on timescales from few months to few 
years (Fig. [3]). This suggests that the frequency of optical 
flares and the ejection rate of jet components are intercon- 
nected. A link between the optical flares and the stationary 
component SI in the radio jet is suggested by a correlation 
between the maxima of optical flares and the characteristic 
epochs tsi of the moving components C4-C12 (Fig. [3]). For 
these nine components, the epochs tsi of separation from the 
stationary feature SI are coincident, within the errors, with 
the maxima in the optical continuum flares. All nine ejec- 
tion events occur within ~ 0.3 yr after a local maximum is 
reached in the intensity of the optical continuum. Consider- 
ing that there are about 14 years of optical data, we estimate 
the probability that maxima of optical flares and ejections 
from SI are associated by chance in the 0.3 yr window. The 
probability of any single ejection event at SI occurring ran- 
domly less than 0.3 yr after a maximum of optical flare is 
about 0.02. The significance of the null hypothesis that all 
nine events happen by chance is w 10 -10 and can be rejected 
at a confidence level of > 99.99 % (for a description of the 
statistical test see Marscher et al. 2002). The significance is 
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still high (~ 10~ 9 ) even without two components C4 and 
C8 for which the gaps in the optical data do not allow the 
timing of the optical maxima although the data is consistent 
with a maximum. 

This is a clear evidence that the passage of moving ra- 
dio knots through the location of the stationary component 
SI and optical continuum flares are physically related. The 
radio events follow the maxima of optical flares with the 
mean time delay of (f» 0.10 ± 0.04) yr. The average distance 
between the region where the optical flare reaches maximum 
and the component SI is measured to be 0.04 pc (or 0.36 pc 
from D) adopting that the knots propagate with the mean 
apparent speed w c during the time period of ~ 0.1 yr. 
Radio events are associated with optical flares occurring on 
time-scales of about a few months to a few years (Fig. [3]). 
The rise of all nine optical flares occurs between two radio 
events, the ejection from component D and separation from 
SI (Fig. [3} suggesting that the generation of variable optical 
emission happens in the innermost part of the jet, between 
D and SI, as a result of propagation of these ejections along 
the jet. The generation of optical flares and its character- 
istics such as the amplitude, timescale and the frequency 
are likely to be related to the properties of the subparsec- 
scale jet: ejection rate, structure and kinematics. Therefore, 
understanding the physical nature of the regions D and SI 
is important for understanding the structure of the central 
engine in 3C 390.3 and location of continuum source(s) pro- 
ducing the variable continuum radiation. 



4 VARIABLE RADIO EMISSION OF THE 

STATIONARY COMPONENTS IN THE JET 

Variations of the flux densities of D and SI components of 
the jet (/d and /si) are shown in Fig. [3] (cyan circles and 
cyan diamonds). The relatively sparse time sampling of the 
VLBI measurements does not allow us to test for a corre- 
lation between radio flares and epochs of ejections. Frac- 
tional var i abilit y amplitude of SI (0.44±0.18), defined as in 
iLu fc Yul | 2001), is about two times higher than that of D 
and the optical emission (0.3 ± 0.17 and 0.24 ± 0.07, respec- 
tively). The optical and radio flux densities of D increase 
by the same factor («3 times) from 1996 to 2008 showing 
a surprisingly similar behaviour over a time period of 14 
years, while the radio flares from SI, varying on scales from 
a few months to a few years, show some correlation with 
the optical flares having comparable timescales. Moreover, 
magnitudes of the optical and radio flares are comparable 
in SI, with the flux density varying by a factor of three for 
optical continuum flux and a factor of four in the radio. It is 
therefore likely that the optical variations on timescales of 
few decades are related to the radio variability of the compo- 
nent D, while the optical flares are coupled with radio flares 
of SI component of the jet. 

We used t he Z-transformed correlation function 
(|Alexanderlll997T l to calculate the correlation between / opt 
and /si as this method can deal with sparsely sampled data. 
The correlation function between / op t and /si have two 
peaks at a timelag around — 50^* days and — 350^ jg days 
(implying that radio flares lead optical flares with time lag 
of about one year), with a correlation coefficient of around 
0.551 

0.22- Correlation function between / op t and /d has a 
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Figure 4. X-ray flux variations in 3C 390.3 and the epoch of jet 
component ejection at D, the base of the jet. Top: Hard X-ray 
fluxes at 2 keV to 20 keV (Gliozzi et al. 2006) and the epoch of 
ejection of C8. Bottom: Soft X-ray fluxes at 0.1 keV to 2 keV 
(Lcighly et al. 1997) and the epochs of ejection of C5 and C6. lc 
error bars are presented for all data. 



flat shape and high correlation coefficient, around 0.85I 07 , 
with time lags from 200 days to 400 days and from -400 
days to -500 days. The present data suggest a correlation 
between variations of optical and radio emission curves but 
uncertainties of the correlation coefficient are large (because 
of the poor radio sampling) for making final conclusions. 
A denser VLBI radio sampling covering a time period of 
several years is needed to improve the correlation observed 
between the radio and optical flares. 



5 PHYSICAL IDENTIFICATION OF THE 
REGIONS D AND SI 

The presence of several apparently stationary features in 
the radio jet of 3C 390.3 may be explained in several dif- 
ferent ways. Accordin g to present models of relativistic jets 
l|Marscher et al.ll2008l . and references therein), the plasma 
material accreted onto central nucleus is collimated into the 
jet. At th e base of the jet , the jet flow either becomes su- 
persomc llDalv et alJll98St i or optically thin (|Konig]| Il98ll : 
lLobanovl 19981) or releases the energy co ntained in the Poynt- 
ing flux ( Romanova fc Lovelace! 1 1996h . The plasma mate- 
rial is then accelerated and collimated into continuous jet 
flow by helical magnetic fields on typical scales of ~ 10 4 
Schwarzschild radii. At the end of the acceleration and col- 
limation zone (< 1 pc), the jet flow may become turbulent 
(due to weakening of the helical magnetic field) and end in 
a stand ing conical shock ass ociated with the radio core of 
the jet (jMarscher et al.ll2008h . 
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The link between the ejection epochs of jet components 
and the dip in the X-ray was found for two radio galaxie s 
3C120 and 3C 111 (|Marscher et al.ll2002l : iMarscherl l2006f ). 
This correlation was interpreted as accretion of the X-ray- 
emitting gas in the disk into the central black holes and 
ejection of a fraction of the infalling matter into the jet 
(Marscher et al.ll2002h . We found a similar correlation be- 
tween the ejection epochs of the components C5 and C8 
and the dip in t he X-ray flux (see Fig. [U) and hardening 
of the spectrum l|Arshakian fc BellonilboOrj ) . The ejection 
of C4 and C8 components happens during the dip in the 
X-ray flux suggesting that the component D is located near 
the black hole. Adopting an anglqj of 50° between the jet 
and the line of the sight, the corresponding de-projected 
distance between D and SI is 0.3 pc/ sin 50° ~ 0.4 pc. 
The average brightness temperature of the component D, 
Tb(D) ~ 1.2 x l O^K, is close to the inverse Co mpton 
limit ~ 5 x 10 10 K (|Kellermann fc Paurinv-Tothlll969l ) while 
the Tb(Sl) ~ 2.2 x 10 10 K of SI component suggests that 
the synchrotron radiation mechanism maintains the approx- 
imate equipartition in energy between the particles and th e 
magnetic field in the region around SI l|Readheadl 1 19941 ) . 
The component D can be identified with the accretion disk 
or hot corona, actual base of the jet or the radio core. 

The high brightness temperature rules out the D be- 
ing the accretion disk or hot corona since the brightness 
temperature for three major accretion disk models is esti- 
mated to be less than 10 10 K. Standard disks are charac- 
terized b y a balance between radia tive cooling and viscous 
heating ( Shakura fc Sunvaevl 1 19731 s ) . For a typical quasar 
with M bh = 1O 8 M and L ~ 10 3S W, the effective tem- 
perature will be T e fl ~ 10 5 K at a radius of r ~ 100r s 
(r s = 2GMbh/c 2 is the Schwarzschild radius for a black 
hole of mass Mbh, where G is the Newtonian gravitational 
constant). However, a hot thin corona (kT e ~ 70 — 250 keV; 
T c ~ 0.8 — 3 x 10 9 K) above the disk layer can be a bright 
source of radiation by means of inverse Compton scatter- 
ing o f UV photons from the disk l|Haardt fc Maraschilll99ll . 
I1993T ). Advection-dominated accretion flow (ADAF) is ra- 
diatively inefficient so t hat most of heat pro duced through 
accretion is preserved l|Naravan fc Yil \l994 ). The ADAF 
model predicts that the mass accretion rate is a factor of 
~ 10~ 4 of the Eddington limit and that t he electron tem- 
perat ure is above 10 9 K for r < 100r s (Manmot o et al.l 
Il997h . The slim disk model, with an extremely high mass 
accretion rate above the Eddington limit, was proposed by 
Hb ramowicz et al.l (|l988T l. The disk is hotter than the stan- 
dard disk because of the high accretion rate. The brightness 
of the blackbody radiation from slim disks is ~ 10 6 K. How- 
ever, as in the case of standard disks, a hot corona may 
produce a radiation with T c ~ 0.8 — 3 x 10 9 K. 

The mean brightness temperature of the D component 
is higher than that of the SI, 7b(D) > T^(S1), suggesting 
that D is the radio core of the jet rather than the base of 
the jet. Were D the core of the counterjet and SI the core 



2 On the assumption that the pattern speed and bulk speed 
of the jet of 3C 390.3 are equal, the jet inclination angle 8 ~ 
50°, bulk Lorentz factor 7 ~ 2 and beaming angle if) PS 
7 — 1 ~ 30° are estimated using the v ariable Doppler factor 
5 = 1.16 (Lahtecnmaki & Valtaoia 1 9991) and the maximum ap- 
parent speed of 1.6 c observed in the compact jet. 



of the approaching jet, we should expect that fo < /si 
because of relativistic (Doppler) beaming effect. The fact 
that fo > /si over the 14 years of monitoring period (see 
Fig. [3]l rules out D being the core of the counterjet. The 
component D should then be associated with the radio core 
of the jet and its emission is non-thermal and generated by 
inverse Compton mechanism (Tt, ~ 1.2 x 10 11 K). Then the 
component SI should be associated with the standing shock 
(a common feature seen in superluminal AGN) downstream 
from the radio core. This is supported by the average Tb ~ 
2.2 x 10 10 K estimated for SI component, which is close to 
the equipartition limit (~5x 10 10 K; Readhead 1994) . 

Alternative scenario assumes that the component D is 
the base of the jet located in the immediate vicinity of the 
central black hole. Then the SI stationary feature should be 
the radio core of the jet. Location of D near the accretion 
disk is supported from the X-ray data (Fig. [4| . Very long- 
term variations (on time-scales of decades) of radio emission 
from D component is not characteristic for the jet emission 
which varies on time scales from months to years. Variations 
from D component on time-scales of decades mimic very 
long-term variations of optical continuum emission (Fig. [3j 
which is probably generated in or above a hot corona at 
a distance ^200 R s above the accretion disk (|Fabianll2004l ; 
iPonti et al.l I2004T ) . The non-thermal radio emission at the 
base of the jet can be generated by synchrotron or inverse 
Compton (7b (D) ~ 1.2 x 10 K) mechanisms by free elec- 
trons of corona accelerated in the rotating and twisted mag- 
netic fields generated by the accretion disk. 

Another possibility is that D and SI are two indepen- 
dent, active nuclei in the binary black hole system. The bi- 
nary black hole scenario cannot explain the correlated, si- 
mult aneojis_jwiahilityo£the red and blue wings of the H/3 
line l|Shapovalova et al]|200ll ). If the two black holes are lo- 
cated in D and SI, the total mass of the binary must ex- 
ceed ~ 5 x 10 9 (P o b s /1000 yr) M . The observed radial ve- 
locity changes indicate a possib le periodicity shorter than 
100 yr (|Shapovalova et aUl200ll ). pushing the total mass of 
the binary to ~ 10 11 Mq, in contradiction with the mass 
of ~ 4 x 1 8 Mo measured u sing the reverberation mapping 
technique dKaspi et al.ll200oT ). The trends observed in the ra- 
dial velocities cannot be reconciled with the orbital motion 
in a binary black hole l|Eracleous et al.| [l997). Thus it is also 
difficult to use a binary system with two black holes located 
in D and SI as an explanation for the observed properties 
of 3C 390.3. 

To understand the radiation mechanism and structure 
of the inner nuclear region of radio-loud galaxy 3C 390.3, it 
is fundamental to identify and locate the D and SI regions. 
This requires multi-frequency VLBI observations to measure 
the spectral index and polarization of stationary components 
of the compact jet. 



6 DISCUSSION AND CONCLUSIONS 

Soldi et al. (20 08) have carried o ut multiwavelength of the 
quasar 3C273 (|Soldi et al.ll2008T ). and have shown that its 
emission is complex due to the presence of many emission 
regions in the central engine. They suggested that variable 
optical emission is non-thermal and it is likely to be gener- 
ated in the base of the jet. Evidence for optical synchrotron 
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emission generated in and near the radio core of the jet 
comes from monitoring of radio a nd optical polarized emis - 
sion of the quasar PKS 0420-014 (jP'Arcangelo et alJliooH ). 
Concurrent changes of the polarization angle at optical and 
radio (43 GHz) indicate that these variable emission is non- 
thermal and originate in the same region around the radio 
core. While the emission of the accretion disk is viewed di- 
rectly for a quasar, it is hidden in radio galaxies by the 
dusty torus. Very weak (or absent) blue bump in 3C 390.3 
|Wamsteker et al1ll997h indicates that the accretion disk is 
obscured and most of variable optical continuum emission 
could be generated in the jet. 

The significant correlation found between the optical 
flares and the epochs tsi can be interpreted in terms of dis- 
turbances (radio knots) propagating from the location of D 
component to the SI stationary component. If the compo- 
nent D is the base of the jet then the following physical 
scenario is likely. A portion of the disk X-ray emitting ma- 
terial accreted onto central nucleus is ejected from the base 
of the jet (in the immediate vicinity from the black hole) at 
the epoch tu- The knot of electrons accelerated to high en- 
ergies radiates a rising emission in X-ray, optical and radio 
as it propagates upstream from the base of the jet. The jet 
becomes visible in radio at tsi because of the e mitting region 
of the jet becoming opticall y thin or denser |Konig]||l98ll : 
iRomanova fc Lovelace! Tl996T l. Alternatively, if the compo- 
nent D is the core of the jet then the rising synchrotron emis- 
sion from X-ray to radio bands can be produced between D 
and SI components by acceleration of the kn ot (shock?) in 
a helical magnetic field IjMarscher et al.ll2008l ) and/or evolu- 
tion of the shock spectrum (or turnover frequency; Valtaoja 
et al. 1992). The minimum Doppler factor D m i n = 1-5 is 
required to increase the optical flux by factor of 2.5 (see 
Fig. |3]) , if the jet flow accelerates from 7 = 1 to 7 = 2 and 
the apparent speed of the jet is 1.6c. This value D m i n = 1.5 
is in reasonably good agreement with the measured v ariable 
Doppler factor ~ 1.2 (|Lahteenmaki fc Valtaoialll999l ). Then 
the component SI should be a standing sho ck formed by 
continuous relativistic flow l| Gomez et al . 1995). The moving 
knot would be compressed by the shock wave thus producing 
variable emission in a wide range of electromagnetic spec- 
tru m as indicated from multi waveband polarization variabil- 
ity l|D'Arcangelo et al]|2007f ). As the disturbance passes the 
standing shock SI, which manifests the end of the accelera- 
tion zone, it expands adiabatically because of low magnetic 
field strength and low electron density, and the optical flux 
drops down on time scales of < 1 yr, while the high fre- 
quency radio emission lasts much longer because low-energy 
electrons have longer lifetimes. In this scenario, the timescale 
and amplitude of non-thermal optical flares depend on kine- 
matics and emitting power of moving perturbations in the 
innermost subparsec-scale region of the jet. 

Variations of the broad-line emission in 3C 390.3 follow 
the optical continuum with some time delay which is esti- 
mated by various a uthors to be in the range from about 
20 days to 100 days dWandel et al.lll999l; iKaspi et al.ll2000l ; 
IShapovalova et al.ll2o"oil ; ISergeev et al.ll2002t ). If the variable 
optical emission is being generated by shocks in the inner- 
most part of the jet then the beamed continuum emission 
from the jet must ionize the gas around the jet thus produc- 
ing the BLR downstream from the 'classical' virialized BLR 
ionized by the optical emission generated in the accretion 



disk. An outflowing, non-virialized BL R can be generated 
in the rotating su brelativistic outflow l|Murrav et al.|[l997l ; 
iProga et al.ll200ol ) surrounding the jet. In BL Lac type ob- 
jects (in which the jet oriented near the line-of-sight) mov- 
ing superluminal components also may produce the short- 
term radio flares, opti cal, X-ray and even gamma-ray flares 
l|Marscher et al.l 120081 ) when passing near the line-of-sight 
of the observer. In the case of 3C 390.3 the jet is inclined at 
« 50°, hence, the Doppler beaming of continuum emission 
is weak towards the observer but the emission is enhanced 
in the direction of the jet and able to ionize the conical out- 
flowing BLR with a half-opening angle 2 ip ~ 7 -1 ~ 30°. In 
radio-loud galaxies there could be two regions of generation 
of variable optical emission, non-thermal emission in the in- 
nermost part of the jet and thermal optical emission in the 
accretion disk. 

The existence of the jet-excited outflowing BLR in 
3C 390.3 will question the as sumption of virilized motion in 
the BLR (|Kaspi et alJl200Ch of all radio-loud AGN, galax- 
ies and quasars, and , hence, the applica bility of the re- 
verberation mapping (|Peterson et al.ll2002h to estimate the 
black hole masses of radio-loud AGN. Time delays and pro- 
file widths measured during periods when the jet emission 
is dominant may not necessarily reflect the Keplerian mo- 
tion in the disk, but rather trace the rotation and out- 
ward motion in an outflow. This can result in large er- 
rors in estimates of black hole masses made from moni- 
toring of the broad-line emission. In the case of 3C 390.3, 
the black hole mass (2.1 x 10 9 Mq) estimated effectively 
from the measure ments near the maximum in the contin- 
uum light curve (IShapovalova et all l200ll ) is significantly 
larger t han the values ((3.5 t o4) x 10 s Mp) rep orted in other 
works (|Wandel et all 1 19991 ; IKaspi et alj|2000h . This differ- 
ence is reconciled by considering the line width (upwhm ~ 
10500 kms^ 1 ) and the time delay (r « 24 days) between 
the optical continuum and line fluxes near the minimum 
of the continuum light curve, which yields Mbh ~ 1-45 x 
10 5 M (cr/lt - day)uf WHM = 3.8 x 10 8 M Q . The possible 
existence of an outflow-like region in a number of radio-loud 
AGN should be taken into account when estimates of the nu- 
clear mass are made from the variability of broad emission 
lines. 

The presence of the jet-excited nonvirial BLR in radio- 
loud AGN is capable of explaining some of the spectral 
characteristics of emission lines. Depending on the orien- 
tation of the jet, the approaching and rotating outflow ma- 
terial in the BLR will imprint prominent signatures on the 
emission lines. At small viewing angles of the jet this BLR 
may produce blue-shifted and single-peaked broad emis- 
sion lines, while non-shifted and double-peak emission lines 
IIEracleous fc Halperni r2003) will be observed at large angles 
of the outflowing BLR to the line of sight. T he narrow emis - 
sion lines will have similar characteristics (|Borosonl l2005h 
being ionized in the approaching subrelativistic outflow by 
the beamed continuum emission of the jet. 

The correlation found for 3C 390.3, between the epochs 
of maxima in optical flares and passing of radio knots 
through SI c omponent, is also confir med for the radio 
galaxy 3C 120 (|Le6n-Tavares et alj|2009h . To check whether 
this correlation is common for other galaxies we conducted 
coordinated long-term radio-optical observations of several 
nearby radio-loud galaxies. A denser VLBI radio sampling 
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covering a time period of several years is required to resolve 
the correlation between variabilities of the radio emission of 
the jet components and the optical continuum, and further 
to constrain the models for the nuclear region of radio-loud 
galaxies. 

The principal results of this work comes from analysis of 
combined radio VLBI (15 GHz), optical/UV and X-ray data 
of the radio-loud galaxy 3C 390.3, and they can be summa- 
rized as follows: 

(i) Structure of the parsec-scale jet: from ten epochs of 
VLBI observations we identified three stationary compo- 
nents (D, SI and S2) and nine moving components (C4-C12) 
of the jet on scales of a few parsecs. Apparent speeds of mov- 
ing components are estimated to be in the range from 0.7 c 
to 1.6 c. 

(ii) We found a new correlation between the local maxima 
in the optical continuum light curve and the epochs at which 
the moving components of the jet pass the stationary radio 
feature SI. Nine radio events follow the peaks of optical 
flares with the mean time delay of 1.2 ±0.4 months. Optical 
flares brighten when radio knots move between stationary 
components D and SI. 

(iii) Identification and location of D and SI regions: anal- 
ysis of available archival X-ray monitoring data and ejection 
epochs of two components, C5 and C8, revealed that the 
ejection of these components from D occurs during the dip 
in the X-ray emission. It is most likely that the component 
D is associated with the VLBI core of the jet (rather than 
the base of the jet) and the component SI is a standing 
shock located at a distance of » 0.4 pc downstream from 
the component D. 

(iv) These results have important implications for the 
structure of the sub-parsec-scale nuclear region of the radio- 
loud galaxy 3C 390.3. We suggest that the relativistic plasma 
of the jet is a dominant source of variable optical continuum 
emission: the bulk of optical continuum optical emission on 
timescales from few months to years are likely to be gener- 
ated by knots of high-energy electrons propagating in the 
innermost 0.4 pc region of the jet between stationary com- 
ponents D and SI. In this scenario, the timescale, ampli- 
tude and frequency of optical synchrotron flares depend on 
energetics, kinematics and rate of ejected radio knots. The 
beamed continuum emission from the jet ionizes a gas in a 
subrelativistic outflow along the jet, which results in a for- 
mation of two non-virialized outflowing BLRs along the jet 
and counterjet. 
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APPENDIX A: PROPERTIES OF THE VLBA 
IMAGES AND MODEL FITS OF 3C 390.3 

Table lATI compares the CLEAN component models (denoted 
"I") and Gaussian model fits (denoted "M") for the VLBA 
data used in the paper. The columns are: Stotai - total flux 
density [mjy/beam]; .Speak _ peak flux density [mjy/beam]; 
Smi-a - minimum flux density [mjy/beam]; x 2 - goodness 
of the fit parameter; er rms . uv - root-mean-square between 
the observed and model visibilities [mjy]. Last row presents 
average ratios between the respective image and model fit 
properties. The ratios are close to unity for the total and 
peak flux densities. The \ 2 parameter and the visibility 
r.m.s. are only slightly larger for the Gaussian model fits, 
which indicates that the model fits represent the structure 
adequately. The higher maximum negative flux density in 
the model fits (column Smin) indicates that that the SNR of 
the Gaussian fits is on average 1.5 times lower than that of 
the VLBI images (due to increased non-Gaussian shapes of 
low-brightness regions). This reduction does not affect the 
fitted values of the component parameters (albeit it does 
increase the parameter errors) as it is related to extended 
emission associated with the underlying flow. 

Properties of individual emitting regions in the radio jet 
are determined by fitting Gaussian patterns to interferomet- 
ric amplitudes and phase closures. Errors of the modelfits are 
estimated by varying the best fit parameters and determin- 
ing 1-ct confidence limits from the resulting \ 2 distribution. 
This method accounts for correlations between parameters 
of adjacent modelfit components and provides conservative 
error estimates. Fidelity of the component sizes and sepa- 
rations obtained from the modelfits is further examined by 
comparing them with resolution limits calculated from the 
SNR of the detections. The resolution limit, 9 m i n , can be 
estimated for a Gaussian modelfit component from 
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where S is the signal-to-noise of component detection and 
&maj, 6min are the major and minor axis of the point-spread 
function derived from the interferometric coverage of the 
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components with sizes 9i t 2 can be fitted adequately to visi- 
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where # m i n ,i and 0min,2 are the minimum resolvable sizes for 
the two components. Table IA2I lists the ratios between the 
measured sizes of and separations of the components D and 
SI and their respective limiting values. For resolved com- 
ponents the ratios between the measured sizes 9 and their 
respective # m i n are larger than unity. For the component sep- 
arations, ratios large than unity imply that the components 
parameters can be reliably estimated from the modelfit. 

The SNR of detection varies from 30 to 900 for SI and 
from 100 to 860 for D and the corresponding resolution lim- 
its are smaller than 0.2 mas for both features. The ratios 
between measured and limiting sizes and separations of D 
and SI indicate that only the size measurements made for 
the first data point for SI and the last one for D are coming 
close to the resolution limits, while all other measurements 
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Table Al. Properties of the VLBA images and model fits of 3C 390.3. 



Epoch ^total ■Speak ^min g rms,uv 

I M I M I M I M I M 



1994.67 


423.9 


422.9 


213.3 


211.1 


—2 


2 


—3.4 


1 


175 


1.186 


253.1 


253.8 


1995.96 


463.1 


466.2 


190.8 


190.8 


— 1 


6 


—2.8 





328 


0.334 


119.2 


120.3 


1996.37 


471.6 


477.6 


192.4 


192.6 


— 1 


3 


—3.0 





453 


0.475 


83.4 


85.5 


1996.82 


389.5 


399.1 


174.2 


174.2 


— 1 


5 


—3.2 


1 


235 


1.504 


42.5 


45.2 


1997.19 


370.0 


377.5 


179.1 


180.4 


— 


7 


—2.2 


1 


533 


1.800 


39.9 


43.0 


1997.66 


418.3 


425.2 


267.4 


267.2 


-1 


8 


-2.2 





575 


0.580 


152.1 


152.5 


1998.21 


392.4 


389.6 


191.2 


191.4 


-1 


1 


-1.3 





429 


0.430 


94.6 


94.7 


1999.85 


313.2 


315.8 


216.9 


218.6 


-1 


1 


-1.3 


1 


103 


1.149 


31.0 


31.4 


2000.65 


310.4 


300.2 


161.8 


167.4 


-5 


3 


-7.3 


1 


234 


1.595 


38.9 


42.2 


2001.17 


270.6 


280.1 


209.4 


207.6 


-1 


2 


-2.6 


1 


250 


1.578 


32.2 


35.4 


2002.13 


301.3 


297.6 


217.0 


217.3 


-0 


9 


-1.8 


1 


094 


1.437 


28.4 


31.6 


2003.02 


268.1 


254.3 


211.0 


212.5 


-3 


3 


-5.3 


1 


275 


1.929 


35.0 


41.2 


2005.13 


299.6 


303.9 


220.3 


217.0 


-3 


7 


-2.8 


1 


013 


0.996 


26.1 


25.5 


2006.11 


344.8 


337.9 


295.9 


287.3 


-1 


7 


-3.1 


1 


074 


1.446 


91.6 


54.7 


2006.31 


340.6 


348.3 


278.8 


279.0 


-3 


1 


-3.4 


1 


127 


1.228 


31.6 


32.5 


2006.54 


390.4 


391.2 


306.5 


305.3 


-3 





-3.5 


1 


049 


1.092 


33.3 


33.7 


2006.75 


379.2 


377.6 


303.0 


303.8 


-3 


2 


-2.8 


1 


052 


1.067 


31.6 


31.9 


2007.01 


405.1 


409.8 


350.5 


350.3 


-1 


6 


-1.8 


1 


109 


1.109 


118.9 


118.8 


2007.63 


444.5 


445.7 


380.6 


380.6 


-1 


6 


-1.8 


1 


131 


1.141 


76.3 


76.3 


2008.41 


468.6 


462.6 


394.3 


394.4 


-2 


1 


-2.7 


1 


189 


1.213 


68.7 


68.5 


2008.60 


429.4 


422.6 


344.0 


341.5 


-1 


9 


-3.1 


1 


138 


1.267 


65.8 


69.2 



< I/M > 0.999 1.000 0.7242 0.906 0.999 



Column designation: Stotal — total flux in the image ("I") and in the respective Gaussian ("M") component models, Speak ~~ peak flux 
density in the respective images and models, S m i n - minimum flux density in the respective images and models, x' 2 ~ goodness-of-fit 
parameter of the respective images and models, cr rms - r.m.s visibility noise for the respective images and models. 



are certainly above these limits, both for the component 
sizes and separations. The limiting ratios and separations 
can be used to estimate the probability of the component 
flux to be correlated: 



/ \ / ^min 



(h 



hi 



\Si-lJ 



This formula gives p c = 1 for two identical Gaussian features 
separated by r m i n . The correlation probabilities (listed in the 
last column of Table IA2[) are smaller than 6 % for all of the 
model fits, implying that the flux densities of the component 
D and SI estimated from these modelfits are uncorrelated. 

This paper has been typeset from a TffjX/ IATgX file prepared 
by the author. 



Table A2. Resolution limits and fidelity of the flux densities in 
the modelfits 



Epoch 


SNR 


8/ #min 


r/r min 






D SI 


D SI 


D-Sl 





1994.67 


433.3 


32.0 


3 


13 





99 


3.8 


0.0540 


1995.96 


454.3 


303.9 


1 


90 


1 


64 


17.6 


0.0094 


1996.37 


388.8 


116.4 


1 


51 


1 


79 


9.1 


0.0182 


1996.82 


481.6 


66.3 


2 


71 


2 


32 


5.3 


0.0608 


1997.19 


861.8 


130.8 


3 


16 


1 


67 


11.0 


0.0156 


1997.66 


475.2 


140.2 


1 


53 


1 


76 


11.3 


0.0115 


1998.21 


535.1 


407.3 


2 


62 


2 


88 


13.5 


0.0401 


1999.85 


568.4 


266.0 


2 


81 


1 


31 


13.4 


0.0205 


2000.65 


115.9 


70.4 


1 


08 


1 


58 


6.1 


0.0431 


2001.17 


508.5 


239.4 


1 


06 


1 


41 


9.7 


0.0127 


2003.02 


632.4 


356.2 


1 


21 


1 


15 


6.6 


0.0272 


2005.13 


406.2 


38.4 


2 


42 





99 


3.8 


0.0428 


2006.11 


499.7 


90.8 


1 


24 


1 


19 


10.6 


0.0046 


2006.31 


455.4 


116.2 


1 


29 


1 


65 


9.3 


0.0122 


2006.54 


681.9 


265.7 


1 


71 


1 


91 


9.4 


0.0249 


2006.75 


573.2 


174.6 


1 


21 


1 


54 


11.0 


0.0088 


2007.01 


688.3 


128.4 


1 


62 


1 


34 


9.0 


0.0097 


2007.63 


1069.1 


351.5 


1 


28 


1 


24 


13.3 


0.0053 


2008.41 


1599.0 


423.5 


1 


89 


1 


02 


10.7 


0.0099 


2008.60 


1633.3 


266.2 


2 


39 


1 


39 


7.0 


0.0245 



Column designation: SNR - signal-to-noise of component detec- 
tion, 0/Bmin - component ratio of the measured and minimum 
resolvable sizes, r/r m ; n - ratio of the measured distances the min- 
imum resolvable separations between the components D and SI, 
p c - probability of the measured flux densities of D and SI being 
correlated. 



